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TRANSLATIONAL MOTION OF BODIES UNDER THE FREE SURFACE OF A 
HEAVY FLUID OF FINITE DEPTH* 


By M. D. Haskind 


In reference 1, entitled "The Two-Dimensional Problem of the 
Vibration of Bodies under the Surface of & Heavy Fluid of Finite Depth," 
the problem was to determine the wave motion of a heavy fluid excited 
by the periodic vibrations of & body of arbitrary shape situated under 
the free surface of the fluid of finite depth; the method of N. E. Kochin 
(reference 2) was used. 


In the present paper, the two-dimensional problem of the wave 
motion produced in a heavy fluid of finite depth by the horizontal 
rectilinear and uniform motion of a solid body of arbitrary shape 
immersed under the surface of the fluid is considered by the same 
method. 


1. Statement of the Problem 


The problem of the translatory motion of & solid body under the 
free surface of a heavy incompressible fluid of finite depth will be 
considered. The case in which the motion of the body occurs with con- 
stant horizontal velocity c will be studied. The motion of the fluid 
will be defined with reference to a moving system of coordinates Oxy 
fixed to the body, the x-axis coinciding with the undisturbed level of 
the fluid and directed along the direction of motion of the body, and 
the y-axis directed vertically upward. 





it will be assumed that the motion of the fluid is potential and. 


E steady relative to the body. From the integral of Lagrange for the 

. pressure within the fluid, 

= : L 

E p- po = pe S pay (1-1) 
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where pg is the atmospheric pressure, p the density of the fluid, 


g the acceleration of gravity; 9(x,y) the potential of the absolute 
motion of the fluid, and v = [ erad 9| the magnitude of the absolute 
velocity of the fluid. 


The function (x,y) is determined from the boundary conditions; 
the flow condition on the wetted contour of the body, 


Se = 6 cos(n,x) on C (1.2) 


where n is the outer normal to the contour C; 
on the free boundary p - Po» and hence 


eS 7 EV - gy a (1.3) 


on the bottom of the channel for y = -hp, the following condition 
applies 


g =0 | (1.4) 


According to the theory of waves of small amplitude, condi- 
tion (1.3) may be linearized. For this purpose the boundary condi- 
tion (1.3) is referred to the x-axis and the term v2/2 neglected. In 
place of condition (1.3), 


$.8yx-o — (1.5) 


It is easily seen that on the free surface the following relation 
holds 


cy(x) = 9+ const (1.6) 
where y is the stream function. In fact, when the stream function of 
the motion of the fluid relative to the body is denoted by Yo, there 
is obtained 

Vo = V - cy 


or 


cy = V - Vo 
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From this relation, equation (1.6) follows ; Since the boundary of 
the fluid in the relative motion is represented by stream lines on which 
Vo is constant. For the free surface, it may be assumed that yo = O. 


Hence, on the free surface, 











| ey(x) = y 
and therefore boundary condition (1,5) assumes the form 
X -vyso for y=0 (1.7) 
where 
v= £ (1.8) 
e 
From condition (1.5) it is seen that the equation of the free sur- 
face will be 
; c [òv 
x) =< 1.9 
j r(x) - 2 È] (1.9) 


2. Fundamental Formulas of the Problem 


The problem may be mathematically formulated as follows. It is 
required to determine the characteristic function w(z) = 9 + iy 


(z = x + iy; i = A/-1), satisfying the conditions: 


l. For O?y7-hgo in the region occupied by the fluid, the deriv- 


ative dw/dz is finite and at infinity for x >+% ; the derivative 
dw/dz vanishes. 


2. On the contour C, the smooth flow condition applies 


E = e cos(n,x) 


5. On the free surface for y = 0, the linearized condition holds 
with regard ‘to the constancy of the pressure 


Re(dw/dz + ivw) = 0 
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. 4. On the bottom of the channel for y = -hg, the following condi-- 


tion holds 
Im dw/dz = 0 


In the region oceupied by the fluid, the point z is taken and 
two contours C4, and C, are drawn, of which Ce contains both the 


point z and the contour C, while the C4, contains the contour C, 


but not the point z (fig. 1). By the formula of Cauchy for a single- 
valued function dw/dz = v(z), 


=a l V(t) ag 1. v(t) à 
v(2) = zu mama ` Zr "RUE (2.1) 
Cy Ca 


where the bar over a letter indicates, as usual, the transition to the 
complex conjugate value. The following notation is introduced 


mG-2h f MEE v=- | IEEE (2 


G 


It is evident that Vj(z) is a holomorphic function in the entire 
plane of the complex variable outside the contour C}, having at infin- 


ity the order zil and capable of being continued analytically in the 
entire part of the complex variable plane which lies outside the con- 
tour C, while Vo(z) is a holomorphic function within the contour Cæ; 


by the extension of which an analytical continuation of this function 
may be obtained over the entire strip O>y>-hp. 


The function Vo(z) may be represented in another form. For this 
it is possible to find a function w(z), which in the strip O>y>-hpo 


has a single pole of the first order € = & + in with residue A/ 21i. 
and which satisfies conditions 1, 5, and 4. 


In fact, for a vortex of strength T, located at the complex point 
€ = £ + in, an expression for the complex velocity was obtained by 
Tikhonov (reference 3) 


eam mr 
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ee eee gene, NU 
orla) = zala -4J ^ Zui(z- t3 Bib) ^ 
Y | i x(n + ha) cos Mec da iho) _ 
E J (eX) ept mo) TT RS o 


" sh Ag(n + bo). 
Y vho - ch? oho 


where Ag. is the real and positive root of the equation 


sin' olz - &+ iho) . (8.5) 


v sh Aho = Ach Aho (2.4) 


For c?« Alea, in all cases where the function to be integrated 


has a singularity, the principal value in the sense of Cauchy is taken 
under the integral. 


For c? > Nao» equation (2.4) has only imaginary roots and the 


fourth term of formula (2.5), which determines the presence of free 
waves, is absent. 


For a source of strength 89 located at the complex point 
b = E inj, the expression of the complex velocity may be obtained in 
the same manner as in the case of a vortex. Without the computations, 
the final result is 


_ Q Q 
mlz) = an(z - €) + 2z(z - 5 + 2iho) s 


ch Aln + ho) sin A(z - E + ihg) 


2 p DYER) eum Mg c aa en a Se E 
O . 


ch Ag(1 + ho) 
vho - ch? Agho 


By the use of expressions (2.3) and (2.5), to obtain the function 
w(z) may be obtained without difficulty. For this purpose, since 
A =T + iQ, the following expression is obtained after simple 
transformations: 


Qv cos Ao(z - E + iho) (2:5) 
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"m | A 
o(2) = zat) tile be Zino) ` 


(Rame ) - A sin A(z - £) 
i (v X) exp(- aho) Tune X ch Aho WA 
O 
sy WA See c (2.6) 


9s vho - ch” Joho 


Here, as in the preceding formulas, the fourth term, which deter- 
mines the presence of free waves, is present only if cĉ< Af&ho- 


When A - v(6£) dë is substituted in the previous formula and 
integration is carried out over the contour (Cj, 


sin A(z - € + 2ihg) 


v(z) = Vy(z ae ors 3i iur. vshAmo- Xen Am; AI 
C1 (0) 


eos Ao(z - € + 2ihg) (z E) a 
i vho - ch? agho SM. x v(t (fen exp(- ano) v im. ch Mg 9 


cos Ag(z - €) : 
my I E (2.7) 


ac- CD 


If both points z and € are situated in the strip O>y>-ho, 
the following equation holds 


S a f in(z - Ẹ + 2iho)| aa (2.8) 
z - Ẹ + 2ihọ A [ | 


With this eguation taken in account, it is found from eguation (2.7) 
that the function Vs(z) can be represented in the form 


ty ys 






EE Tinea 


TEE 


e tees 





mc 
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Va(z) = - af ep f - i exp EG -Ẹ + 2ino)] + 
d UC Ci i 0 


sin Mz - € + 2ihg) cos Ao(z - € + 21ng) | _ 
(v +A), exp(- aho) Wah Rig = heh Kay T rate d i: 


EN * js s sin Mz - £) oa uc ROME LOS KEER | 
af +2) exp(- Aho) V sh ho - X ch Ahg a o- w vho - chê Agho ab 


- (8.9) 


The conjugate complex functions are introduced for real X 


HQ) = v(t) exp - ixtat, HO) JA v(6) exp ixtadt 


Cy Cy (2.10) 


By an interchange in equation (2.9) of the order of integration, 
and by simple transformations ; there is readily obtained 


Vo(z) EM A) ex[ a + 2ino)] + 


(v X) exp(- Aho) 


Z(v sh Xhg - X ch Ahg) (xc- A) exp [56 * 2ino)| - EQ) ex[ - in(z + 21m) z 
H(A) exp idz + H(- A) exp(- 12) dA - 


niv 


Z(vh - ch® Aghg) (sc Ao) exp ido (z+ 2iho) + E(Ag) exp [- iolz + zimo) | s 


H(1o) exp idoz - H(- Xo) exp(- 092) (2.11) 


It is of interest to find the character of the waves that remain 
behind the moving body. For this purpose the asymptotic expression of 


. the. complex. velocity is first obtained for x>- e in the case of a 


vortex and source. In reference 5, the asymptotic expression of the 
complex velocity in the case of a vortex is of the form 
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ha i PEE 
(pl, = - 2TV INT sinAg(z - & + ihg) (2.12) 


In a similar manner, the asymptotic expression of the complex 
velocity is obtained in the case of a source. Without the computations, 


the final result is 


| " ch Ao(n + ho) 


(ag) a p - yng = ene AOG cos ho(z -č + iho) (2.13) 


For the function o(z), having a polarity with residue A/axi ; the 
following asymptotic expression is obtained: 
A cos A(z - € + 2ihg) - A cos Xo(z - t) 


2 (2.14) 
v hg - ch Xogho 


(5) a cd - - iv 


Setting A= v(t) a¢ and integrating over the contour C, yields 
the asymptotic expression of the function v(z) = dw/dz: 


z) iv [i : ; 
a = o OR H(- x ) exp iX (z + 2ih ) + 
(3 X -> 2(vhg x ch? zobo) 9 i 9 0 


H(A) exp [ - Dolz + 21h9)] - H(Ag) exp Doz - R(- Xo) exp(- Doz)| 
(2.15) 


Finally, from the formula 


c dw 
x) = < Re[{ = 
y( ) 8 (E) 


it is readily found that for x--œ sinusoidal waves of length 2x/A9 


are formed behind the amplitude of which, after some simple transforma- 
tions, may be represented in the form 


a WA a H(A) exp àgbo - H(- Xo) exp - Xoho | 


(2.16) 








e. 


IP tL E 


moet cec 





AA RL cc M re OT aaa 
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3. Formulas for Determining the Forces 


The forces acting on the contour C are now computed. The lift 
force of the contour is denoted by P, the resistance by R, and the 
moment of the forces on the contour about the origin by M. These 
forces will be computed by the formulas of Chaplygin-Blasius: 


WA -a fn de ane g [sot az 
Co. ZEE (e l 


2 2 l (3.1) 


where Cp is an arbitrary contour, situated in the region 0> y> -ho 
and containing the contour C; and vo(z) is the complex velocity in 


the relative motion obtained by superposing on the absolute flow a 
uniform motion of the fluid with velocity c in the direction of the 
negative x-axis. Thus, 


vo(z) - V4 (z) + Vo(z) - Cc 
where the contour Cy is chosen to lie between C and Co. 


Formulas (3.1) do not take into account the buoyancy force of 
Archimedes, equal to gpS, and it moment, equal to -gpSx,, where S 


is the area that bounds the contour C, and x, is the abscissa of the 
center of gravity of this area. 


The following integral is now computed: 


J de dz dA dz [eem - c)? dz + dA - c) dz 
C2 C2 C2 C2 


But the first and second integrals on the right are equal to zero 
because the function V4(z) is holomorphic outside the contour Co 


and has at infinity a zero of at least the first orđer, while the func- 
tion V2(z) is holomorphic within the contour Cs. Hence, 


J= 2 [wat dz - z f'on + Vo) dz = VA dz - e [sl dz 
C C C C 


2 2 2 
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The velocity circulation about any contour that contains the con- 
tour C is denoted by T, so that 


r= v(z) dz 


therefore 


P-iR=- [mn dz + pcr X2) 
C2 


By the use of expressions (2.2) and aly, the following expression is 
obtained 


V4 (z) Vo(z) dz = af Da à "MI A) exp in(z + 2ihg) + 
Co 


MA H(- A) exp iA(z + 21hg) - H(A) exe[ - id(z + 21n9)] 
z- expl- kho) ¥ sh Aho - Ach Abo : 





y + 3 H(A) exp iÀàz - H(- A) exp(- =) . 
V expt kho) v sh kio - ACH Aho x 


aay E(- Ag) exp Dofz + 2iho) + H(X\o) exp [ - ixo(z. + 21ho]] " 


2 vho - ch" agho 
xiv H(Ag) exp idoz + E(- Ao) exp(- iñoz) "m 
2 vho - chê Agbo 


Since the point 0, which belongs to the contour Cj, lies within 
the contour “Cs, with an interchange in the order of integration and 


by the following formula, 


1 etià z az = eHAE 
ani ze 
C2 


There is obtained 
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» Vi(z) Vg(z) dz = al BA scu exp(- 23hg) + 
Co 0 


at- ») |? exp (- 23hg) - | H(A) | 2: exp Ato 
v sh Abo - Ach Aho 


y 





À exp(- X ho) 


xiv B(- Ao)l? exp(- oho) + | HQ) |? exp Zoho - 2H(Ag) B(- X o) 
2 vho - ch? Agho 


Hence, formula (2.3) assumes the form 


3.3) 


P - iR = pcr - $- j |I8C2)]? exp(- 2Xn9) + 
0 


2391? ES - |E)? dh, 
(v +2) exp(- Aho) EM ee aoa m D o] a 





< vho - ch? Agho (5.4) 


Separating the real and imaginary parts and adding to P the 
Archimedes force, not taken into account by the Chaplygin-Blasius for- 
mula, results in 


1 | P = per - & f [x a)? expl- 23ho) + (v + X) exp(- Aho) x 
A i 


1 Iat- XI* esc ano) E 4, m[sog Cao). 5 uu 


i | 2(v sh 3ho - Ach dbo 2(Vhg - ch? zobo) 





gi JO)! exp cho + | 8(- 39)]? exp(- Agbo) - 2 Ref HO) &C- Ao)) 
Aa vho - ch? Aghg 


n 
pi = 
|| 


(3.6) 
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Formula (5.6) may be given another form, namely 
hus pv. RO) exp Agbo - H(- Xo) exp(- Roho)! 4 
4 ch? X obo - Vho 


It can be readily shown that the total resistance of the underwater 
wing consists only of the wave resistance. In fact, by the following 
well-known formula for computing the wave resistance in the case of & 
fluid of finite depth, 


(5.7) 


2X ^h, 
Ql 2 -00 
R = 7 pga (a zu xL) (3.8) 


„anā with the value of the amplitude a from formula (2.16), for- 
mula (3.7) is obtained after some transformations. 


The moment of the acting forces on the contour C is now com- 
puted. When the moment of the Archimedes force is taken into account, 


M = - gpSxe + Re EC + Ve(z) - cl? az (3.9) 
CL 


This expression is computed in an entirely similar manner to the 
computation of the expression P - iR. 


For very large absolute values of z the following expansion can 
be employed 


vi) = za CS. “ain | TO de. ute 
el 


C 








and, hence, 


Further, 





OUR 
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| * l I z(V3 Š De dz 20 
peres | "T | n 5 


and therefore, 
M = - gpSx, + wf Vi (Vp - c) dz 
. l ; Co 


or, since the function Vo(z) is holomorphic within the contour Cg 


M = - gpSx, - pc se [ate dz + p AA vı (z) Vo(z) az 
C Cs 


2 (3.10) 


It is noted that 


ueta: YA EF) exp(- ang) at 
C2 


The integrals in formula (3.10) are computed in the same manner as in 
the expression (3.3), and as a result there is obtained the formula 


= - goSx, - pe Re (int (0j] +O SEYA [e (- 2%) El- X) exp(- 2X9ho) + 
O 


(v + X) expl- A ho) 


Z(v sh Xhg - A sh Aho) (maha X) expl- 2Xho) + BU) 30) exp 2Aho - 
z'(- A) BO) = H(- A) m) | an - 


l T a (zc Ao) E(- Ao) exp(- 2Aoho) - H'(Ao) H(Ro) exp 2dgho - 


H'(- Ao) H(A) + B(- Ao) sos) ` (5.11) 
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Formulas (3.5), (3.7), and (3.11) in the limiting. case for hoe 
agree with the formulas obtained by Kochin in reference 2. 


The function H(A) in formulas (3.5), (3.7), and (3.11) does not 
depend on the contour C4, and for example, the contour C or some 


other contour which contains the contour C may be taken for the con- 
tour of integration. Moreover, the value of the function H(A) does 
not change if, instead of the complex velocity v(z) of the absolute 
motion, the complex velocity of the relative motion vg(z) is taken, 


because these two functions differ by a constant c. The properties 
of the function H(A) will be used in the following section. 


4. Examples 


In the preceding sections expressions were found in terms of the 
function H(A) of a number of important magnitudes, namely, the ampli- 
tude of the waves formed, the wave resistance, the lift force, and the 
moment of the forces acting on the contour. Thus, the function 


H(A) [ve exp(- 1X6) at gi exp(- ix£) (4.1) 
C l C 


plays a fundamental part for the problem under consideration. In order 
to compute this function, it is necessary to know the expression for 
the complex velocity, i.e., the solution of the hydrodynamic problem. 
In ease the relative depth of the submerged contour C is sufficiently 
large, however, & good approximation is obtained if, in place of the 
function v(z), there is substituted in formula (4.1) the expression 

of the complex velocity which corresponds to the motion of the con- 
tour C in an infinite fluid. 


Several examples of such an approximate solution of the problem 
will be considered 


1. The motion of a circular cylinder. - The circular cylinder of 
radius b, situated at the depth h under the free surface of the 
fluid, is assumed to move with constant horizontal forward velocity c, 
since the circulation about the contour of the cylinder has & given 
value IT. In this case, the characteristic function for the infinite : 
fluid is known: 


cb? r ; 
wW (z) ae gee Inte ih) 


Gcr CM E MET eer e mine 


eS 


aA EE 


€: 


"EE 




















Ti cine: 


mereatur 
€ 


Ex 





fomes 


Ó— 


IIS caer aa 
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Hence, 


v(z) = [PE + tile) i (4.2) 


By formula (4.1) the function H(A) is now constructed: 


cb? 


z r“ " ] ^ E 
H(A) = (ee mye + seen] exp - ilz dz 
C 


Since the contour C contains one singular point z = -ih, there 
is obtained by the theorem on residues 


H(A) = (T + 2xcb?1) exp - Ah (4.3) 


With the use of formula (3.7), the expression for the wave resist- 
ance of the cylinder is obtained 


"m mo sh Xo(ho - h) + 2xe*9b? ch Ag(hg - x)? ea 


ch? Aghg - vho 


and by the use of formula (5.5) the expression for the lift force of 
the cylinder is obtained 


P = pc - pre ds pceb?T " pc? b 


4x(ho - B) 2(no - n)? ^ 2(ng - n) 


i J WA emus (T + 4x?c?p 3?) sh 2A(ho - h) + 4xcb?T^ ch 2A(hg - n) 
vc exp(- 
o 


9s veh Aho - À ch Aho dd + goS = (4.5) 


The integral component of this formula may be computed by the 
method of mechanical quadratures. In the limiting cases v=O and 
v =, this component can be very accurately computed. Moreover, if 
this integral component is considered as a function of the parameter 
a = l/(vho) = c?/(ghg), it can be shown that for a = 1 this component 


suffers a discontinuity. In the particular case when the radius b of 
the cylinder is taken equal to zero, i.e., when the motion of a vortex 
under & free surface is.considered, formulas (4.4) and (4.5) lead to 
the expressions established by Tikhonov. It is noted further that for- 
mulas (4.4) and (4.5) have been derived on the assumption that c@< gho. 


For c?» gho, no free waves are formed behind the cylinder and the wave 
resistance R is equal to zero. 
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For the moment of the forces exerted by.the fluid on the cylinder 7 
the following expression is obtained by formula (3.11): 


H'(- Ao) H(- Ao) exp(- ohio) - E! (ko) H(A) exp Doh | 
vho - ch? Roho 


ov E' CC Ao) EQ) - B(- Ro) E' (ào) 
4 vho - ch? Agho 





.. BY 
Meo, 





But from equation (4.5), it is evident that 
H'(A) = - hE(3) + 2xcb? exp(- Ah) 
H'(- X) = hH(- 3) - 2xcb? exp Ah 
Hence, after simple transformations, 


P sh? Ag(hg - h) + seb^Ao sh 2Ao(ho - h) 
ch? Aghg - vho 


(4.6) 


M = hR - 2xpcb?v 


The point of intersection with the y-axis of the resultant force 
on the body is determined by the formula 





M 2xcb? 
yo- R =- h+ 2 (4.7) 
T + 2zxcb^ÀAo cth Ao(ho - h) 
It is evident that for R>O this resultant never passes through 
the center of the cylinder. 


2. Motion of an elliptic cylinder. - An ellipse, having & center 
at the depth h and having axes Za and 2B directed parallel to 
the axes of coordinates x and y, is allowed to move with & constant 
velocity c in the direction of the x-axis. The circulation T is, 
for simplicity, taken equal to zero. In this case, the flow of an 
infinite fluid about the contour C is determined with the aid of an 
auxiliary variable and the formula 


1 2 
z=- ih+3 Na? - p? EBE ve alee (nee) 
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where r= Ala. + B)/(a - B) andlul =r is the equation of the cir- 


cle in the u-plane which corresponds to the contour of the ellipse C. 
The exterior of this circle corresponds to the exterior of the ellipse. 
The following function is set up: 


H(A) a eat ikz) dw = 
WC 


L(g) WA me Rcg) 


ful =r 


When the substitution u = iv is made, there is obtained 


2 
_ H(A) =- E Na? - BÊ exp(- Ah) Jt + =) exp x Ala? - pê (v - 3) dv 


Ivi zr 


But by the theory. of Bessel functions it is known that 





l 
a A & 2 xp 3 ( E i) = Jaz) 
Jv] =r 


hence, , 
H(A) = ze Ala? - p2 ex (9) | 2. GIG = 87) 4 s (WNGE ~ Be) } 
From the formula 

J.3(z) = - Iz (z) 


and the value of r, the following expression is obtained 


H(A) = ancp AJ SE 2L exp(- àh) ONZ - p2) (4.8) 


The computation is restricted to the wave resistance. By 
formula (3.7), 
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sado deg NB. d. pec aay | 
R= a eee” CB ae RE Jji^(oNa^ - B^) | (4.9) 


From this formula, it follows that for certain Ag and, therefore, 
for & certain velocity c <Negho, the wave resistance is equal to zero; 


i.e., the amplitude of the waves formed behind the moving.body becomes 
zero. This will be tbe case if the following relation is satisfied: 


Ao Na? - p? = Sk (k21,2, . . .) 


where s, is the positive root of the Bessel function Jj(s). The 
first root of this function is 


$81 = 3.832 


Since the parameter V = g/ c? is connected with Xo by the 
equation 2 


c? 


the first velocity at which the wave resistance becomes zero is deter- 
mined by the formula 


5.832hg 


e = 0.51 AJgANa? - B2? th ——— (4.10) 
Aa? - p? 
Moreover, 
5.832 
aE Rede. < 1 
Als - Be 
hence, 


e « 0.51 Nena? - g2 (4.11) 


In a similar manner a number of other examples may be considered. 
Moreover, as in reference 2, it is possible in this case to set up a 
functional equation for determining the function H(A) and the values 
of the circulation [ from the condition of the finite velocity at the 
sharp edge. These equations may be obtained by the same method. Their 
final form will be somewhat more complicated as compared with the case 
of the infinite fluid. 
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Translated by S. Reiss x ao X& poseen 


National Advisory Committee 
for Aeronautics 
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